Systems of reaction-diffusion equations with 
spatially distributed hysteresis 

Pavel Gurevich, Sergey Tikhomirov 



1. Setting of the problem. Reaction-diffusion equations with spatially dis- 
tributed hysteresis were first introduced in [I] to describe the growth of a colony of 
bacteria (Salmonella typhimurium) and explain emerging spatial patterns of the bac- 
teria density. In [U [2], numerical analysis of the problem has been carried out, however 
without rigorous justification. First analytical results were obtained in [31 H] (see 
also [5j), where existence of solutions for multi- valued hysteresis was proved. Formal 
asymptotic expansions of solutions were recently obtained for some special case in [5]. 
Questions about the uniqueness of solutions and their continuous dependence on initial 
data as well as a thorough analysis of pattern formation remained open. In this paper, 
we formulate sufficient conditions that guarantee existence, uniqueness, and continu- 
ous dependence of solutions on initial data for systems of reaction-diffusion equations 
with discontinuous spatially distributed hysteresis. Analogous conditions for scalar 
equations have been considered by the authors in [71 E] - 

Denote Q T = (0, 1) x (0,T), where T > 0. Let U C R k and V C R l (fc, I G N) be 
closed sets. We assume throughout that (x,t) G Q T , u(x,t) G U, v(x,t) G V. 

We consider the system of reaction-diffusion equations 

u t = DAu + f(u,v,W(Z ,u)), 

(1) 

v t = g(u } v,W(£ ,u)) 
with the initial and boundary conditions 

u\t= Q = (fi{x) 1 v\ t=0 = il)(x), (2) 
u x \x=o = u x \x=i = 0. (3) 

Here D is a positive-definite diagonal matrix; W is a hysteresis operator which maps 
an initial configuration function ^o(^) (G {1,-1}) and an input function u(x, •) to an 
output function W(£o(x), u(x, As a function of (x,t), W(^q,u) takes values in a 

set W C M m (m G N). Now we shall define this operator in detail. 

Let r^r^ C U be two disjoint smooth manifolds of codimension one without 
boundary (hysteresis "thresholds"). For simplicity, we assume that they are given by 



Ja(u) = and 7/?(«) = with Vj a (u) 7^ and V^p(u) ^ 0, respectively, where 7 a and 
7/3 are C°°-smooth functions (in the general situation, atlases can be used). 

Denote M a = {u G IA : 7 a («) > 0}, = {u G U : 7 /3 (w) < 0}, M a/3 = {u e U : 
7 a (u) < 0, 7/3 (w) > 0}. Assume that M a nT f3 = and Afy D T Q = (Fig. [Lj). 

Next, we introduce locally Holder continuous functions (hysteresis "branches") 

Wi : D(Wt) = M a U M Q/3 ->> W, VF_! : D(W_!) = M p U M Q/3 W. 

We fix T > and denote by C r [0,T) the space of functions which are continuous 
on the right in [0,T). For any ( G {1, —1} (initial configuration) and u G C([0,T]; U) 
(input), we introduce the configuration function 

C : {1, -1} x C([0,T]-U) -+ C r [0,T), ((t) = C(Co,«o)(*) 

as follows. Let X t = {s G (0,t] : %(s) G T Q U T^}. Then C(0) = 1 if w (0) G M a , 
C(0) = -1 if «o(0) G M/3, C(0) = Co if «o(0) G M a/3 ; for t G (0,T], = C(0) tiX t = 0, 
((t) = \ ii X t ^ and wo(maxXi) G T a , ((t) = — 1 if X t ^ and Mo(maxX t ) G 

(Fig. ED. 




Figure 1: Regions of different behavior of hysteresis W 

Now we introduce the hysteresis operator W : {1,-1} x C([0,T];U) —> C r [0,T) 
by the following rule (cf. P, [TO]). For any initial configuration Co G {1,-1} and input 
wo G C([0,T]-U), the function W(( ,u ) : [0,T] -)• W (output) is given by 

W r (Co > «o)(*) = Wc (t) (u (t)) > (4) 

where ((t) is the configuration function defined above. 

Now assume that the initial configuration and the input function depend on spatial 
variable x G [0, 1]. Denote them by £o(x) and u(x,t), where u(x, ■) G C([0,T];U). Let 
u(x, 0) = (f(x) for x G [0, 1]. 



Condition 1. The initial data <p(x) and ^o(x) are consistent in the following sense. 
For any x G [0,1], we have ^o(x) = 1 if <p(x) G M a , ^o(x) = — 1 if (p(x) G Mp, and 

£ (x) e {1,-1} if <p(x)eM aP . 

Now, using ((H), we can define the function 

W^oOn), «(*,-))(*), (5) 

where x is treated as a parameter. We call this function spatially distributed hysteresis 
and the function £(x,t) = ((£o(x),u(x, -))(t) its spatial configuration. 

Definition 1. We say that a function (p G C' 1 ([0, 1];U) is transverse with respect to 
a spatial configuration £o(^) if it is consistent with £o(x) and the following holds. If 
(p(x) G T a (r^) and the vector (f'(x) is tangent to r a (Tp) at some point x G [0, 1], 
then £o(X) = 1 (~1) m a neighborhood of x. 

Condition 2. The initial function (f(x) is transverse with respect to £o(x). 

Remark 1. Suppose that ip(x) G T a . Then the tangency of the vector (p'(x) to the 
manifold r a can be analytically written as 



0. 



where (•, ■) is the inner product in 



In this paper, we shall deal will spatially distributed hysteresis whose spatial con- 
figuration has finitely many discontinuity points in x for each t. 

Condition 3. There exist points = b < ~&i < ■ ■ ■ < ~&m < &m+i = 1 (M > 1) such 
that £ (x) = {-iy {or i-lY +1 ) for x G ft, b l+1 ), i = 0, . . . , M. 

Now we consider time-dependent functions u(x } t) such that u,u x G C(Q T ;U). 

Definition 2. We say that a function u is transverse on [0, T] (with respect to a spatial 
configuration if, for every fixed t G [0,T], the function u(-,t) is transverse with 

respect to the spatial configuration £(■,£). 

Definition 3. We say that a function u preserves spatial topology (of a spatial config- 
uration t)) on [0, T] if there is M > such that, for t G [0, T], there are continuous 
functions = b (t) < bi(t) < • • • < b M {t) < b M +i{t) = 1 such that £(x,t) = (— 1) % (or 
respectively) for x G (&<(*), h+i(t)), i = 0, . . . , M. 



The functions bi(t) play a role of free boundary, which has much in common with 
free boundary in parabolic obstacle problems (see, e.g., [TTT [T2"] and the references 
therein). However, in our case, the behaviour of bi(i) is defined differently. 



2. Functional spaces for solutions. We denote by L q (0, 1), q > 1, the standard 
Lebesgue space and by M^(0, 1) with natural I the standard Sobolev space. For a 
noninteger I > 0, denote by VF'(0, 1) the Sobolev space with the norm 



w « \Jo Jo |x-y|i+9('-ffl) y y ' 
where [/] is the integer part of I. Further, we introduce the anisotropic Sobolev spaces 
W^{Qt) with the norm ( \\u(-, t)||^ 2 dt + / Q T <ZtV * and W^HQt) with 

the norm |M|loo(Qt) + H^tlUocCQr)- We denote by C X (Q T ), < A < 1, the Holder space. 

For the vector-valued functions, we use the following notation. For example, if 
u(x,t) G U and each component of u belongs to W^ ,1 (Qt), then we write u G 
W 2 /{Q T] U). 

Throughout the paper, we fix q and A such that q G (3, oo) and A G (0, 1 — 3/g). In 
particular, this implies that u,^ G C X (Q T ;U) whenever u G Wg ,1 (QT'M) (see Lemma 
3.3 in [13, Chap. 2]). 

To define the space of initial data, we use the fact that if u G Wjj' 1 (Qt'M) , then 
the trace u\ t=Q is well defined and belongs to Wq 2//<? ((0, 1);U) (see Lemma 2.4 in [T3"l 
Chap. 2]). Moreover, one can define the space W^ ((0,1);U) as the subspace of 
functions from Wq 2 ^ g ((0, 1);W) with the zero Neumann boundary conditions. 

We assume that tp G W^~ N 2/q ((0, 1); U) and ip G ^((0, 1); V) in @. 

Definition 4. A pair (u(x,t),v(x,t)) is called a solution of problem (JTJ) — (J3J) (in Qt) if 
it G Wq ,1 (Q T ;U), v G M / ^ 1 (Qr;V); W(£ ,u) is measurable with respect to (x,t); u and 
f satisfy equations ([I]) for a.e. (x,t) G Qt! the first condition in (J2I and conditions (J3J) 
are satisfied in the sense of traces; the second condition in (j2J) holds in the following 
sense: v(-,t) ->■ ifj(-) as t ->■ in Loo((0, 1); V). 

The solution is called transverse (preserving spatial topology) if u(x } t) is transverse 
(preserving spatial topology). 

It follows from this definition that any solution u(x,t) belongs to C(Q T ;U). There- 
fore, the function W(£o, u) is well defined by (03) and belongs to L oq (Qt; VV). 



3. Assumptions on the right-hand side. First, we assume the following. 

Condition 4. The functions f(u,v,w) and g(u,v,w) are locally Lipschitz continuous 
%nU x V x W. 

Next, we formulate dissipativity conditions for / and g. 

In the following condition, we denote by U^, [i > 0, parallelepipeds in U with the 
edges parallel to respective coordinate axes such that (f(x) G £Y M for all x G [0, 1]. 

Condition 5. There is a parallelepiped U and, for each sufficiently small (i > 0, there 
is a parallelepiped and a locally Lipschitz continuous function f^{u,v) such that 

1. \f^(u,v)\ converges to uniformly on compact sets in U x V as fj, — >■ 0, 

2. At each point u G OUq D D(W±i), v G V, the vector f(u,v,W±i(u)) + f^{u,v) 
points strictly inside Uq. 

3. At each point u G dU^ D D(W±\), v G V, the vector f(u,v,W±i(u lx )) + f^{u,v) 
points strictly inside for all G U^. 

To formulate the assumption on g, we fix some numbers T > and M > 0. 
Consider arbitrary functions ip(x) G V, Uq(x, t) G £Y, Wq(x, t) G W such that 

II^IUoo((o,i) ; v) < M , \\uq\\ Loo{Qto . u) < M , \\w \\ Loc{Qtq;W) < M Q . 

It was proved in [T41 Theorem 1, p. Ill] that the system of ordinary differential 
equations (with x being a parameter) 

Vt = g(u (x,i),v, w (x,t)), v\ t=0 = ip(x) (6) 

either has a unique solution v(x,t), \\v\\loo{Qt ;V) < 00 > or there is a maximal existence 
time T max G (0,T ) such that \\v(-, t) \\ Loo ((o,iy,v) -> oo as t ->• T max . 

Condition 6. For any T > and M > 0, t/iere zs a compact Vq = Vq(T , M ) C V 
which does not depend on ift,Uo,Wo such that the solution v(x,t) of system (ED satisfies 
v(x,t) G Vo for all (x,t) G Qt - 

Remark 2. In particular, the uniform boundedness of v holds if \g(u, v,w)\ < 
w)\v\ + B(u,w), where A(w, w) and B(u,w) are bounded on compact sets (see 
Example [ID. However, if V 7^ R, one must additionally check that i> does not reach 
the boundary of V in finite time. To fulfill Condition El one could alternatively assume 
the existence of invariant parallelipipeds for g (similarly to Condition [5]) . 



Example 1. The assumptions of the present paper hold for a model describing a 
growth of a colony of bacteria (Salmonella typhimurium) on a petri plate (see [H 
[2]). Let Ui(x,t) and U2(x,t) denote the concentrations of diffusing buffer (pH level) 
and histidine (nutrient), respectively, while v(x 1 t) denote the density of nondiffusing 
bacteria. These three unknown functions satisfy the following equations: 



supplemented by initial and no-flux (Neumann) boundary conditions. In those equa- 
tions, Di, D 2 , a, <2i, d2 > are given constants, while W(^q,u) is the hysteresis op- 
erator. In this example, we have U = {u £ M 2 : Ui,u 2 > 0}, V = [0, oo), 
W = [0, oo). The hysteresis thresholds r a and Tp are the curves on the plane given by 
7«(«) := Ui — a a /u2 — b a = and jp(u) := U\ — ap/u2 — bp = 0, respectively, where 
tt a ,ci/3,b a ,b/3 > are some constants (Fig. [[]); the hysteresis "branches" are Wi(u) = 1 
and W-xiu) = 0. 

4. Main results. In this section, we formulate the main results of the paper. 
They can be proved by the methods developed in [71 [S]. 
In what follows we assume that Conditions [TH6] hold. 

Theorem 1 (local existence). There is a number T > (depending on ip and ip) such 
that the following holds for problem (jTJ) — (ED in Qt- 

1. There is at least one solution. 

2. Any solution is transverse and preserves spatial topology. 

Theorem 2 (continuation). Let (u,v) be a transverse topology preserving solution of 
problem ([T]) -([3D in Qt for some T > 0. Then it can be continued to an interval 
[0,T max ), where T max > T has the following properties. 1. For any t < T max , the 
pair (u,v) is a transverse solution of problem (CO) -([3D in Qt - 2. Either T max = oo, 
or T max < oo and (u,v) is a solution in QT ma x> but u(-,T max ) is not transverse with 
respect to £{-,T max ). 





Note that a solution need not be topology preserving inside its maximal interval of 
transverse existence. However, if a transverse solution preserves spatial topology and 
is unique, then it continuously depends on the initial data (p, and £o- 

Now we discuss the uniqueness of solutions. Let Uq be the set from Condition [51 

Condition 7. There are numbers K > and a G [0, 1) such that 

K 

\Wi{u) - Wi(w)| < - — rr~ ' ~, — tttt-|w - u\ Vu, u G M a n M a/3 , 

(7/3 Wr + (7/3(w)) CT 

X 

\W-x{u) -W- X {u)\ < \u-u\ Vu,wG Mp Pi M a p. 

Theorem 3 (uniqueness). Assume additionally that Condition^ holds. Let (u,v) and 
(u,v) be two transverse solutions of problem ([I])-(jHj) in Qt for some T > 0. T/jen 
(w, v) = (u, v). 
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